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We observe photons and neutrinos from stars. Based on these observations, complemented by 
measurements of cosmic rays energies and composition, we have been able to constrain several models 
for the Big Bang and for stellar evolution. But that is not enough. We also need to help this effort 
with laboratory experiments. We are still far from being able to reproduce stellar environments in 
a terrestrial laboratory. But in many cases we can obtain accurate nuclear reaction rates needed for 
modeling primordial nucleosynthesis and hydrostatic burning in stars. The relevant reactions are 
difficult to measure directly in the laboratory at the small astrophysical energies. In recent years 
indirect reaction methods have been developed and applied to extract low-energy astrophysical 
S-factors. These methods require a combination of new experimental techniques and theoretical 
efforts, which are the subject of this short review. 



I. ASTROPHYSICS: WHAT WE CAN AND 
WHAT WE CAN'T DO 

A. Hot plasmas on Earth 

Evidently, we cannot reproduce in the laboratory con- 
ditions existing during the Big Bang and during stellar 
evolution. But efforts to reproduce such conditions on a 
limited scale on Earth are underway. A good example 
are experiments being carried out at the National Igni- 
tion Facility (NIF) in Livermore. In this facility the in- 
tense energy of 192 giant laser beams is focused on a small 
spherical pellet containing a few milligrams of fusion fuel, 
typically a mix of deuterium and tritium. The energy 
heats the surface of the pellet into a plasma, exploding 
off its surface, driving the remaining portion of the target 
inwards, and compressing it into a high density. A shock 
wave travels towards the center of the compressed fuel 
from all sides, further heating and compressing it so that 
fusion reactions will occur and release energy, creating 
temperatures and pressures similar to those that exist 
only in the cores of stars and giant planets and inside 
nuclear weapons [T]. 

Another example is ITER, a large-scale international 
laboratory located in France that aims to demonstrate 
that it is possible to produce commercial energy from fu- 
sion. ITER is based on the "tokamak" concept of mag- 
netic confinement, in which the plasma is contained in a 
doughnut-shaped vacuum vessel. A mixture of deuterium 
and tritium is heated to temperatures of 150 million °C, 
forming a hot plasma. Strong magnetic fields are used to 
keep the plasma away from the walls. From 50 MW of in- 
put power, the ITER machine is designed to produce 500 
MW of fusion power. ITER runs on a predicted 15 billion 
euros building cost, whereas NIF already costs roughly 
US$ 5 billions. So, these are not cheap machines at all. 
It is very hard to reproduce conditions within stars. And 
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the prospects of generating energy for commercial use 
with similar projects in the future are still uncertain. As 
for helping us understanding features of the Big Bang 
and of stellar evolution, ITER will not be able to tell 
us much. It will mainly access questions on atomic and 
material science associated with confining a plasma at 
huge temperatures within a vessel and the interactions 
of the plasma with the walls of the vessel. While it is 
undeniable that this experiment will fill a knowledge gap 
needed for further developments in science, it will not 
answer crucial questions of relevance for astrophysics [2] . 
Maybe ITER-2 will, if we can afford it. 

To avoid sounding too negative, I mention that NIF 
has a good plan to provide results on atomic and nuclear 
physics for stellar evolution. I give a couple of exam- 
ples. In the theoretical modeling of stellar evolution one 
relies strongly on calculations of radiation propagation 
through hot stellar plasmas. The coefficient entering the 
radiation propagation equation is called the "opacity". 
It accounts for the interaction of photons with atoms 
and effects such as excitation and ionization of ground- 
state, excited, or ionized atomic species present in the 
medium. For many years we have relied on a huge effort 
to calculate all of the atomic physics needed for stellar 
evolution codes in the form of opacity tables [Sj. Stel- 
lar modelers have not questioned much the reliability of 
such tables, as one simply can't do better than that. But 
it would be a great knowledge improvement if we could 
effectively "measure" opacity in the laboratory. The NIF 
X-ray opacity platform will enable detailed studies of the 
radiative properties of hot dense matter over a photon 
energy range of 200 - 10,000 eV, also important in as- 
trophysics |3]. It will allow benchmarking opacities used 
in the standard solar model and in stellar equilibrium 
codes (relevant to exoplanet habitability assessment) and 
absorption/emission spectroscopy of photoionized plas- 
mas scaled to black hole and neutron star accretion-disk 
conditions. The development of pulsed power and high 
power lasers opens a brand new perspective for the study 
of opacities in several dense plasmas including modeling 
of the atmospheres of very cool white dwarf stars [5]. 
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II. NUCLEAR REACTIONS 

Stars are powered by nuclear reactions at very low en- 
ergies and, in many situations, at very high densities. 
Usually, one needs to know what happens during binary 
encounters between nuclei (a counter-example is the cel- 
ebrated triple-a reaction). The effects of the environ- 
ment electrons are still a disputed research topic. But 
the main problem here is really to know the reaction 
rates at the energies required for stellar modeling. For 
example, in our Sun the reaction ^Be(p, 7)^B plays a ma- 
jor role for the production of high energy neutrinos from 
the /3-decay of ^B. These neutrinos come directly from 
the center of the Sun and are ideal probes of the sun's 
structure. John Bahcall frequently said that this was the 
most important reaction in nuclear astrophysics [6| . Our 
knowledge about this reaction has improved considerably 
due to new radioactive beam facilities. Another example, 
the reaction ^^C(a,7)^^0, is extremely relevant for the 
fate of massive stars. It determines if the remnant of a 
supernova explosion becomes a black-hole or a neutron 
star. These two reactions are just two examples of a 
large number of reactions which are not yet known with 
the required accuracy needed in astrophysics. 

NIF has reported the first cross section and spec- 
tral measurements of the T(t,2n)^He reaction that is an 
important mirror reaction to the '^He(^He,2p)'*He reac- 
tion (which is part of the proton-proton chain in hydro- 
gen burning stars). These direct measurements, which 
were conducted at energies inaccessible by conventional 
accelerator-based techniques, are not affected by electron 
screening. Measurements of the differential cross section 
for the elastic n-'^H and n-^H scattering at 14.1 MeV have 
also been published [7|. The accurate determination of 
this reaction rate is essential for understanding how the 
fuel is assembled in an implosion, and for the demon- 
stration of thermonuclear ignition and net energy gain 
at NIF. It also opens the door for planning the use of 
NIF and other laser powered facilities to obtain informa- 
tion on nuclear reaction rates at the energies occurring 
in stars. 

The extremely low cross sections for reactions induced 
by charged particles and the inherent difficulty to obtain- 
ing reaction cross sections induced by low energy neu- 
trons leads to enormous hurdles to develop reliable stel- 
lar evolution models and computer codes. Chains of low 
energy nuclear reactions lead to complicated phenomena 
such as nucleosynthesis, supernovae explosions, and en- 
ergy production in stars. An example is that approxi- 
mately half of all stable nuclei observed in nature in the 
heavy element region, A > 60, are produced during the 
"r-process" . The exact site of the r-process is not known, 
but one believes that it occurs in environments with large 
neutron densities leading to neutron capture times much 
smaller than the beta-decay half-lives, Tn ^ T/j, of the 
nuclei involved. The most neutron-rich isotopes along 
the r-process path have lifetimes of less than one sec- 
ond; typically 10~^ to 10~^ s. Cross sections for most of 



the nuclei involved are hard to measure experimentally. 
Sometimes, theoretical calculations of the capture cross 
sections and of the beta-decay half-lives are the only 
source of input for r-process modeling. 

Nucleosynthesis in stars is also complicated by the 
presence of electrons. They screen the nuclear charges, 
therefore increasing the fusion probability by reducing 
the Coulomb repulsion. Evidently, the fusion cross sec- 
tions measured in the laboratory have to be corrected 
by the electron screening when used in a stellar model. 
This is a purely theoretical problem as one can not ex- 
actly reproduce the conditions at stellar interiors in the 
laboratory. At least for now. 

A simpler screening mechanism occurs in laboratory 
experiments due to the bound atomic electrons in the 
nuclear targets. This case has been studied in great de- 
tail experimentally, as one can control different charge 
states of the projectile-l-target system in the laboratory 
[51-fT^. The experimental findings disagree systemati- 
cally by a factor of two or more with theory. This is 
surprising as the theory for atomic screening in the lab- 
oratory relies on our basic knowledge of atomic physics. 
At very low energies one can use the simple adiabatic 
model in which the atomic electrons rapidly adjust their 
orbits to the relative motion between the nuclei prior to 
the fusion process. Energy conservation requires that 
the larger electronic binding (due to a larger charge of 
the combined system) leads to an increase of the rela- 
tive motion between the nuclei, thus increasing the fusion 
cross section. As a matter of fact, this enhancement has 
been observed experimentally. The measured values are 
however not compatible with the adiabatic estimate [S]- 
I12j . Dynamical calculations have been performed, but 
they obviously cannot explain the discrepancy as they 
include atomic excitations and ionizations which reduce 
the energy available for fusion. Other small effects, like 
vacuum polarization, atomic and nuclear polarizabilities, 
relativistic effects, etc., have also been considered |13) . 
But the discrepancy between experiment and theory re- 
mains [lo t fTa ] . 

A possible solution of the laboratory screening prob- 
lem was proposed |14[ I15j . Experimentalists often use 
the extrapolation of stopping power tables [15] to obtain 
the average value of the projectile energy due to stopping 
in the target material. The stopping is due to ionization, 
electron-exchange, and other atomic mechanisms. How- 
ever, the extrapolation is challenged by theoretical calcu- 
lations which predict a lower stopping. Smaller stopping 
was indeed verified experimentally TIT . At very low ener- 
gies, it is thought that the stopping mechanism is mainly 
due to electron exchange between projectile and target. 
This has been studied in Ref. 17] in the simplest possible 
situation: proton-t-hydrogen collisions. The calculated 
stopping power was added to the nuclear stopping power 
mechanism, i.e. to the energy loss by the Coulomb re- 
pulsion between the nuclei. The obtained stopping power 
is proportional to v", where v is the projectile velocity 
and a = 1.35. The extrapolations from stopping power 
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tables predict a smaller value of a. Although this result 
seems to indicate the stopping mechanism as a possible 
reason for the laboratory screening problem, the theoret- 
ical calculations tend to disagree on the power of v at low 
energy collisions [18] , 

Another calculation of the stopping power in atomic 
He++He collisions using the two-center molecular orbital 
basis was reported in Ref. [IQj . The agreement with the 
data from Ref. [18] at low energies is excellent. The 
agreement with the data disappears if nuclear recoil is 
included. In fact, the unexpected "disappearance" of the 
nuclear recoil was also observed in Ref. [ID]. This seems 
to violate a basic principle of nature, as the nuclear re- 
coil is due to Coulomb repulsion between projectile and 
target atoms [THj. After several attempts, sometimes 
with elaborate theoretical models, little theory activity 
in this field has been reported. Some models have been 
praised as solving the stellar screening problem (see, e.g. 

I beheve that this praise is more due to the use of 
quantum-field theoretical tools, which tends to impress 
low-energy experimentalists and theorists who know lit- 
tle about those theoretical techniques. The fact is that 
the present situation on screening of nuclear reactions is 
confusing. Either experimentalists are publishing wrong 
analysis, or all aspects of theory might not have been 
considered yet [55]. 

I have discussed above a few ongoing efforts in nu- 
clear astrophysics, some expensive investments to get so- 
lutions, as well as some problems that might not even be 
possible to study on Earth. Actually, much of the knowl- 
edge required for understanding the physics of the Big 
Bang and of stellar evolution can be accessed by means 
of indirect methods in nuclear physics. The goal of these 
theoretical methods, and laboratories that use them, is 
to take a detour following a much harder work of putting 
pieces together from several, sometimes seemingly unre- 
lated, experiments. In this article I will review some of 
these methods and what has been and can be accom- 
plished with them. 



III. UNDERSTANDING FUSION CROSS 
SECTIONS 

A. Fusion 

All approaches to understand fusion reactions in- 
volve two prongs: a) Calculate an ion-ion (usually one- 
dimensional) phenomenological potential (Wood-Saxon, 
proximity, folding, Bass, etc.) using frozen densities, or 
microscopic, macroscopic-microscopic methods using col- 
lective variables (CHF, ATDHF, empirical methods), and 
b) employ quantum mechanical tunneling methods for 
the reduced one-body problem (WKB, IWBC), incorpo- 
rating quantum mechanical processes by hand, including 
neutron transfer and excitations of the entrance chan- 
nel nuclei (CC). Only for very light ions, involving nuclei 
lighter than oxygen it is possible to devise more micro- 



scopic methods, based on binary nucleon-nucleon inter- 
actions, to obtain the fusion reaction cross sections of 
interest for nuclear astrophysics [55] . 

Fusion cross sections can be calculated from the equa- 
tion 

gf{E) = i,\^Y.^U^\)P,(E), (1) 

where E is the center of mass energy, A ^Jh^ /2mE IS 
the reduced wavelength and £ = 0, 1, 2, • • • . The cross 
section is proportional to ttA^, the area of the quantum 
wave. Each part of the wave corresponds to different im- 
pact parameters having different probabilities for fusion. 
As the impact parameter increases, so does the angular 
momentum, hence the reason for the 2£-\-l term. Pi{E) 
is the probability that fusion occurs at a given impact pa- 
rameter, or angular momentum. The barrier penetration 
method (BPM) assumes that fusion occurs when the par- 
ticle (with mass m) penetrates the Coulomb barrier and 
Pp is calculated in a one-dimensional potential model, 
e.g. by using the WKB approximation or alike. From 
ai = 7rA^(2£ -I- l)Pi one can calculate the average value 
of t from {t{E)) — 'Y^i^o't/ '}2iCri and many other rele- 
vant quantities. Sometimes, for a better visualization, or 
for extrapolation to low energies, one uses the concept of 
astrophysical S-factor, redefining the cross section as 

apiE) = ^S{E) exp [-27rr,{E)] , (2) 

where rj{E) = ZiZ2e^ /hv, with v being the relative ve- 
locity. The exponential function is an approximation to 
Po{E) for a square- well nuclear potential plus Coulomb 
potential, whereas the factor \/E is proportional to the 
area appearing in Eq. [l] 

In order to use Eq. ^ one needs the nucleus-nucleus 
potential. This is a badly known beast. It includes 
the effects of non-fusion channels, which might be hardly 
known. As it cannot be calculated from first principles, 
one adds an imaginary part to the real potential and it 
becomes much more than a beast; something really ab- 
noxious. Some have tried to tame this thing from first 
principles. But, except for few heroic attempts, we seem 
to have given up. We just fit whatever we can fit and we 
get whatever parameters of a potential function we can. 
Then we simply call it the "optical potential" . 

B. Many reaction channels 

The situation is worse, as Eq. [T]does not work in most 
situations. A good example is shown in figure [l] (top), 
taken from Ref. [245 • Only by including coupling to other 
channels, the fusion cross sections can be reproduced. In 
coupled channels schemes one expands the total wave- 
function for the system as 

= ^a,{a,qk)'f){a,qk), (3) 
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FIG. 1: Top - Fusion cross section of ''"Ni+'^^Ni as a function 
of the center of mass energy 24; . The dashed (sohd) curve 
is a BPM (coupled-channel) calculation. Bottom - Transfer 
probabilities for multineutron transfer in ^''Zr+'^^Ca [44] , 



where (j) form the channel basis, a is a dynamical vari- 
able (e.g., the distance between the nuclei), and qk are 
intrinsic coordinates. Inserting this expansion in the 
Schrodinger equation yields a set of CC equations in the 
form 



dak 
da 



(4) 



where U is whatever potential couples the channels k 
and j and Ea = Ea — Ea is some sort of transition 
energy, or transition momentum. In the presence of con- 
tinuum states, continuum-continuum coupling (relevant 
for breakup channels) can be included by discretizing the 
continuum. This goes by the name of Continuum Dis- 
cretized Coupled-Channels (CDCC) calculations. There 
are several variations of CC equations, e.g., a set of dif- 



ferential equations for the wavefunctions, instead of using 
basis amplitudes. Coupled channels calculations with a 
large number of channels in continuum couplings, is one 
of the least controllable calculations. Anything can hap- 
pen because of the phases of matrix elements: the cou- 
plings can add destructively or constructively, depending 
on the system and on the nuclear model. Such suppres- 
sions or enhancements are difhcult to understand. 



C. Radiative capture 

For reactions involving light nuclei, only a few channels 
are of relevance. In his CclS6, cL real potential is enough for 
the treatment of fusion. For example, radiative capture 
cross sections of the type n + x a + j and ttL (tt — 
E, {M) =electric (magnetic) L-pole) transitions can be 
calculated from (see, e.g., ^5] ) 



d.c. 



= const. X \ IIOttlII lbjb)\ , 



(5) 



where Ot^l is an EM operator, and {Icjc Ibjb) is a 

multipole matrix element involving bound (&) and contin- 
uum (c) wavefunctons. For electric multipole transitions 

/■OO 

{lcjc\\OEL\\lbjb) ^ const. X / dr r^Ub{r)uc{r), (6) 

Jo 

where Ui are radial wavefunctions. The total direct cap- 
ture cross section is obtained by adding all multipolarities 
and final spins of the bound state {E = Enx), 



■iE) = J2iSF)j,at-XiE) , 

L,Jb 



(7) 



where {SF)j^ are spectroscopic factors. 



D. Asymptotic normalization coefficients 

In a microscopic approach, instead of single-particle 
wavefunctions one often makes use of overlap integrals, 
/f,(r), and a many-body wavefunction for the relative mo- 
tion, Uc{r). Both Ib{r) and Uc{r) might be very compli- 
cated to calculate, depending on how elaborated the mi- 
croscopic model is. The variable r is the relative coordi- 
nate between the nucleon and the nucleus x, with all the 
intrinsic coordinates of the nucleons in x being integrated 
out. The direct capture cross sections are obtained from 
the calculation of af^ oc | (/f,(r)||r^yL||*c(r)) p. 

The imprints of many-body effects will eventually dis- 
appear at large distances between the nucleon and the nu- 
cleus. One thus expects that the overlap function asymp- 
totically matches (r — )■ oo), 

hyf) = (^1 for protons, 

'2k 



Ib{r) = 6*2^ —A';, +i/2(Kr) for neutrons, (8) 
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FIG. 2: Top - Potential model calculation for the reac- 
tion ^Li (p, 7)*Be. Experimental data are from Ref. [26j . 
Bottom - Single-particle model calculations for the reaction 
^Be(p, 7)*B. The dashed-dotted line is the calculation for the 
Ml resonance at Ecm ~ 0.63 MeV and the dotted line is for 
the non-resonant capture. Experimental data are from Refs. 
The total S factor is shown as a solid line. 



where the binding energy of the n + x system is related to 
K by means of Ef, = fi^K^ /2mnx, Wp^q is the Whittaker 
function and is the modified Bessel function. In Eq. 
[sj Ci is the asymptotic normalization coefficient (ANC). 

In the calculation of cr£- above, one often meets the 
situation in which only the asymptotic part of Ib{f) and 
^ c{t) contributes significantly to the integral over r. In 
these situations, Uc{r) is also well described by a simple 
two-body scattering wave (e.g. Coulomb waves). There- 
fore the radial integration in cr^- can be done accurately 
and the only remaining information from the many-body 
physics at short-distances is contained in the asymptotic 
normalization coefficient Ci, i.e. cr'i^^j^ oc C^. We thus 
run into an effective theory for radiative capture cross 
sections, in which the constants Ci carry all the informa- 
tion about the short-distance physics, where the many- 
body aspects are relevant. It is worthwhile to mention 
that these arguments are reasonable for proton capture 
at very low energies, because of the Coulomb barrier. 

As the overlap integral, Eq. |8j asymptotically becomes 
a Whittaker function, so does the single particle bound- 
state wavefunction . If we call the single particle ANC 
by 6i, then the relation between the ANC obtained from 



experiment, or a microscopic model, with the single par- 
ticle ANC is given by {SF)ibf = Cf. This becomes clear 
from Eq. [7] The values of {SF)i and bi obtained with 
the simple potential model are useful telltales of the com- 
plex short-range many-body physics of radiative capture 
reactions [25]. 

Many reactions of interest for nuclear astrophysics in- 
volve nuclei close to the dripline. To describe these re- 
actions, a knowledge of the structure in the continuum 
is a crucial feature. One basic theoretical problem is to 
what extent we know the form of the effective interac- 
tions for threshold states. It is also hopeless that these 
methods can be accurate in describing high-lying states 
in the continuum. In particular, it is not worthwhile to 
pursue this approach to describe direct nuclear reactions. 



E. Resonating group method 

One immediate goal can be achieved in the coming 
years by using the Resonating Group Method (RGM) or 
the Generator Coordinate Method (GCM). These are a 
set of coupled integro-differential equations of the form 

E / rfV[i?l^(r,r')-i?iV„^i^(r,r')]<?a'(r')=0, (9) 

a' 

where H^^,{r,r') = {"if Aia,r)\H\-ii B{a' ,r')) and 
N^^,{r,r') = (*A(a,r)|*B(Q!',r')). In these equations 
H is the Hamiltonian for the system of two nuclei (A 
and B) with the energy E, '^a,b is the wavefunction 
of nucleus A (and B), and ^^(r) is a function to be 
found by numerical solution of Eq. [9] which describes 
the relative motion of A and B in channel a. Full an- 
tisymmetrization between nucleons of A and B are im- 
plicit. Modern nuclear shell-model calculations, includ- 
ing the No-Core-Shell-Model (NCSM) are able to pro- 
vide the wavefunctions 'i>A,B for light nuclei [23]. But 
the Hamiltonian involves an effective interaction in the 
continuum between the clusters A and B. Overlap inte- 
grals of the type lAa{r) = {'^'A-al'^A) for bound states 
has been calculated within the NCSM. This is one of 
the inputs necessary to calculate S-factors for radiative 
capture, Sa ^ \{ga\OEM\lAa)\'^ , whcre Oem is a corre- 
sponding electromagnetic operator. The left-hand side 
of this equation is to be obtained by solving Eq. [9] For 
some cases, in particular for the p-|-^Be reaction, the dis- 
tortion caused by the microscopic structure of the cluster 
does not seem to be crucial to obtain the wavefunction 
in the continuum. The wavefunction is often obtained 
by means of a potential model. The NCSM overlap in- 
tegrals, lAa, can also be corrected to reproduce the right 
asymptotics [34l[35], given by lAa{r) cx VF_^,;+i/2(2fcor), 
where rj is the Sommerfeld parameter, / the angular mo- 
mentum, fco = ^2iJ,Eo/h with /i the reduced mass and 
Eq the separation energy. 

A step in the direction of reconciling structure and re- 
actions for the practical purpose of obtaining astrophys- 
ical S-factors, along the lines described in the previous 
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FIG. 3: World data on ^Be(p,7)*B compared to theoretical 
calculations. 



paragraph, was obtained in Ref . [34j . The wavefunctions 
obtained in this way were shown to reproduce very well 
the momentum distributions in knockout reactions of the 



type ^B+A 



'Be+X. The astrophysical S-factor for 



the reaction ''Be(p,7)®B was also calculated and excel- 
lent agreement was found with the experimental data in 
both direct and indirect measurements [SH [35] . The low- 
and high-energy slopes of the S-factor obtained with the 
NCSM is well described by the fit 



SniE) = (22.109 eV.b) 



1 -I- 5.30£: + 1.65£;2 + 0.857£;3 



i + E/o.mb 



(10) 



where E is the relative energy (in MeV) of p-|-^Be in their 
center-of-mass. This equation corresponds to a Fade ap- 
proximant of the S-factor. A subthreshold pole due to the 
binding energy of is responsible for the denominator 
ESI 137] • Figure [| show the world data on ■^Be(p,7)**B 
compared to a few of the theoretical calculations. The 
recent compilation published in Ref. [22] recommends 
5i7 = 20.8 ± 0.7 (expt) ± 1.4 (theor) eV b. 



IV. DIRECT REACTIONS AND THE ROLE OF 
RADIOACTIVE BEAMS 

A. Transfer reactions 

Transfer reactions A{a, b)B are effective when a mo- 
mentum matching exists between the transferred particle 
and the internal particles in the nucleus. Thus, beam en- 
ergies should be in the range of a few 10 MeV per nucleon. 
Low energy reactions of astrophysical interest can be ex- 
tracted directly from breakup reactions A+a — > b+c+B 
by means of the Trojan Horse method (THM) |38j. If 
the Fermi momentum of the particle x inside a = {b + x) 
compensates for the initial projectile velocity Va, the low 



energy reaction A + x = B + c is induced at very low 
(even vanishing) relative energy between A and x. To 
show this, one writes the DWBA cross section for the 
breakup reaction as 



dflhdflrdEh 



where 



y^?]m(ka,kb,kc)5'faY";m(kc) 



--<xi x\Xa 'Pbx > 



The threshold behavior Ex for the breakup cross sec- 
tion aA+x^B+c = {'n:/kl)Y,i{2l + l)\Si:c? is well 
known: since \Six\ ^ exp(— 27r77), then oa+x^b+c ~ 
(1/fc^) exp(— 27r?7). In addition to the threshold behav- 
ior of Six, the breakup cross section is also governed 
by the threshold behavior of fi{r), which for r — > oo 
is given by //^ (A:^^)^/^ exp(7r?7) K2i+ii0, where 
Ki denotes the Bessel function of the second kind of 
imaginary argument. The quantity ^ is independent 
of kx and is given by ^ = {8r /asY^^ , where as = 

/mZAZxe^ is the Bohr length. From this one ob- 
tains that {d^ /dQbdflcdEb){Ex — )■ 0) « const.. The co- 
incidence cross section tends to a constant which will in 
general be different from zero. This is in striking con- 
trast to the threshold behavior of the two particle reac- 
tion A + x = B + c. The strong barrier penetration effect 
on the charged particle reaction cross section is canceled 
completely by the behavior of the factor Ti„i for 77 — ^ oo. 
Basically, this technique extends the method of trans- 
fer reactions to continuum states, very successful results 
using this technique have been reported [39, 40^. 

Another transfer method, coined as Asymptotic Nor- 
malization Coefficient (ANC) technique pTl - H5] relies on 
fact that the amplitude for the radiative capture cross 
section b + x — > a + 7 is given by 



M 4",(rb^)|0(rbx)|^.^+^(rbx) >, 



where 



^bx 



'< 4>a{ib, £.x, TCx,y,)\(l)x{^x)4>b{£.b) > 



is the integration over the internal coordinates ^b, and 
^x, of 6 and x, respectively. For low energies, the 
overlap integral I^^ is dominated by contributions from 
large rhx- Thus, what matters for the calculation of 
the matrix element M is the asymptotic value of I^x ^ 
Clx W^va,i/2{'2i^bxrbx)/rbx, where Cj}^ is the ANC and W 
is the Whittaker function. This coefficient is the prod- 
uct of the spectroscopic factor and a normalization con- 
stant which depends on the details of the wave function 
in the interior part of the potential. Thus, C^^ is the 
only unknown factor needed to calculate the direct cap- 
ture cross section. These normalization coefficients can 
be found from: 1) analysis of classical nuclear reactions 
such as elastic scattering [by extrapolation of the exper- 
imental scattering phase shifts to the bound state pole 



7 



in the energy plane], or 2) peripheral transfer reactions 
whose amplitudes contain the same overlap function as 
the amplitude of the corresponding astrophysical radia- 
tive capture cross section. 

To illustrate this technique, let us consider the proton 
transfer reaction A{a, b)B, where a — b + p, B = A + p. 
Using the asymptotic form of the overlap integral the 
DWBA cross section is given by 



Jsja 



{CI,? 



Ap 



/3f 



hp 



where a is the reduced cross section not depending on 
the nuclear structure, /3f,p [Pap) are the asymptotic nor- 
malization of the shell model bound state proton wave 
functions in nucleus a{B) which are related to the cor- 
responding ANC's of the overlap function as {C^^? = 
^bpPbp- Here S^p is the spectroscopic factor. Suppose 
the reaction A{a, b)B is peripheral. Then each of the 
bound state wave functions entering a can be approxi- 
mated by its asymptotic form and a cx P^pl^hp- Hence 
da/dn = E,SC%p?iCSp?RBa where Rsa = a/Pl^f^l 
is independent of /3^p and (3^^. Thus for surface reac- 
tions the DWBA cross section is actually parameterized 
in terms of the product of the square of the ANC's of the 
initial and the final nuclei (C'^p)^(C^p)^ rather than spec- 
troscopic factors. This effectively removes the sensitivity 
in the extracted parameters to the internal structure of 
the nucleus. One of the many advantages of using trans- 
fer reaction techniques over direct measurements is to 
avoid the treatment of the screening problem [321 US] ■ 

But do we really understand transfer reactions well 
enough? Let us take an example from literature [44] . 
Assuming that a in Eq. [4] is simply the time t, and using 
the first-Born approximation (i.e, taking ak ^ aoSko), the 
amplitude to excite the channel (t>k from an initial channel 
(j)o is given by = -ih J (^o|C^l0fc) exp[i{Ek - Eo)t/h]. 
The Born approximation can be applied to transfer reac- 
tions. The probability to transfer a nucleon in nucleus A 
from channel a to a nucleon in nucleus B in channel f3 is 
given by 



P, 



Pa 



ih / dtFpa (R-) exp 



. Er — En)t 



^ -)-(•••) 



(11) 

where R is the nucleus- nucleus distance and F^ct(R) is 
the from factor given by 

F^^iR) = J dVe*Q"-0/3(R + r) (Vi - {U)) (/)„(r), (12) 

where Q is the momentum transfer in the reaction, U 
is the total (optical) potential, and Vi is the potential 
of the nucleon with one of the nuclei. Why not V2? In 
the literature, using Vi (V2) goes by the name "prior" 
("post)-form. It has been shown in the past that the 
post and prior forms of breakup and transfer reactions 
lead to the same result. 



In figure [T] (bottom) one sees the probabilities for 
multi-nucleon transfer in ^^Zr-|-'"'Ca, as a function of 
the closest approach distance D — (ZiZ2e^/2_E)[l -|- 
l/sin(0/2)]. Transfer is most likely to occur when 
the nuclei are at their closest point, D. The tunnel- 
ing probability depends exponentially on this distance, 
Ptr/ sin(0/2) exj>(—2aD). This approximation arises 
from Eqs. [TT] and [T2j If one neglects correlations, two- 
nucleon transfer probabilities are given in terms one- 
nucleon transfer probabilities: P2n = {Pin?- For three- 
nucleon transfer P^n = PinP2m and so on. These are 
shown by the straight lines in figure [1] (Right). All seems 
to work well, except that one needs an enhancement of 
a factor 3 to get P2n from theory [45]. That is what 
happens when theorists do not know what to do [44] . 



B. Intermediate energy Coulomb excitation 

At low-energies, the theory of Coulomb excitation is 
very well understood [46|. A large number of small cor- 
rections are now well known in the theory and are nec- 
essary in order to analyze experiments on multiple ex- 
citation and reorientation effects. At the other end, the 
Coulomb excitation of relativistic heavy ions is character- 
ized by straight-line trajectories with impact parameter 
b larger than the sum of the radii of the two colliding 
nuclei [48] . It was also shown that a quantum theory for 
relativistic Coulomb excitation leads to modifications of 
the semiclassical results |4£|. In Refs. t5(ij51j the inclu- 
sion of relativistic effects in semiclassical and quantum 
formulations of Coulomb excitation was fully clarified. 

Recently, the importance of relativistic effects in 
Coulomb excitation of a projectile by a target with charge 
Z2, followed by gamma-decay, in nuclear reactions at in- 
termediate energies was studied in details. The Coulomb 
excitation cross section is given by 



dcTi. 



■f 



da 

dn 



167r^Z|e^ 



ttA/j, 



B{TTX,h^If 

(2A + 1)3 



|^(^A,m)P, (13) 



where B(TrX,Ii — >■ //) is the reduced transition proba- 
bility of the projectile nucleus, ttA = El, E2, Ml, ... is 
the multipolarity of the excitation, and /i = —A, —A + 
1,...,A. 

The relativistic corrections to the Rutherford formula 
for {d(j/d^l)^^ has been investigated in Ref. [17]. It was 
shown that the scattering angle increases by up to 6% 
when relativistic corrections are included in nuclear col- 
lisions at 100 MeV/nucleon. The effect on the elastic 
scattering cross section is even more drastic: up to 13% 
for center-of-mass scattering angles around 0-4 degrees. 

The orbital integrals S{'KX,fi) contain the information 
about relativistic corrections. Inclusion of absorption ef- 
fects in S{TrX,fi) due to the imaginary part of an opti- 
cal nucleus-nucleus potential where worked out in Ref. 
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[50] , These orbital integrals depend on the Lorentz fac- 
tor 7 = (1 — w^/c^)~^/^, with c being the speed of light, 
on the multipolarity ttA/i, and on the adiabacity param- 
eter ^(6) = Ufib/"fv < 1, where ujfi = {Ef — Ei) /h is 
the excitation energy (in units of h) and b is the impact 
parameter. 

Ref. [S2] has shown that at 10 MeV/nucleon the rel- 
ativistic corrections are important only at the level of 
1%. At 500 MeV/nucleon, the correct treatment of the 
recoil corrections is relevant on the level of 1%. Thus the 
non-relativistic treatment of Coulomb excitation |46| can 
be safely used for energies below about 10 MeV/nucleon 
and the relativistic treatment with a straight-line tra- 
jectory |48] is adequate above about 500 MeV/nucleon. 
However at energies around 50 to 100 MeV/nucleon, ac- 
celerator energies common to most radioactive beam fa- 
cilities, it is very important to use a correct treatment 
of recoil and relativistic effects, both kinematically and 
dynamically. At these energies, the corrections can add 
up to 50%. These effects were also shown in Ref. [5T| for 
the case of excitation of giant resonances in collisions at 
intermediate energies. 

A reliable extraction of useful nuclear properties, like 
the electromagnetic response (B(E2)-values, 7-ray angu- 
lar distribution, etc.) from Coulomb excitation experi- 
ments at intermediate energies requires a proper treat- 
ment of special relativity [5H [S3]. The dynamical rel- 
ativistic effects have often been neglected in the analy- 
sis of experiments elsewhere (see, e.g. [55]). The effect 
is highly non-linear, i.e. a 10% increase in the veloc- 
ity might lead to a 50% increase (or decrease) of certain 
physical observables. A general review of the importance 
of the relativistic dynamical effects in intermediate en- 
ergy collisions has been presented in Ref. [511 [SS] ■ 



C. The Coulomb dissociation method 

The Coulomb dissociation method is quite simple. The 
(differential, or angle integrated) Coulomb breakup cross 
section for a -I- ^ — > b + c + A follows from Eq. [13] It 
can be rewritten as 



Eq. 14 is based on first-order perturbation theory. It 



b+c 



(14) 



where lo is the energy transferred from the relative mo- 
tion to the breakup, and a^^,^ b+d^) photo nu- 
clear cross section for the multipolarity ttA and photon 
energy lu. The function F'^^ depends on lu, the relative 
motion energy, nuclear charges and radii, and the scatter- 
ing angle — {0,(1)). F'^^ can be reliably calculated [15] 
for each multipolarity ttA. Time reversal allows one to 
deduce the radiative capture cross section b + c — > a -I- 7 
from cr^^a _j. j_|_^(a;). This method was proposed in Ref. 
[57] and has been tested successfully in a number of re- 
actions of interest for astrophysics. The most celebrated 
case is the reaction ''Be(p,7)*B [SB], followed by numer- 
ous experiments in the last decade (see e.g. Ref. [55]). 



also assumes that the nuclear contribution to the breakup 
is small, or that it can be separated under certain ex- 
perimental conditions. The contribution of the nuclear 
breakup has been examined by several authors (see, e.g. 
[60]). has a small proton separation energy (« 140 
keV). For such loosely-bound systems it had been shown 
that multiple-step, or higher-order effects, are impor- 
tant [ST]. These effects occur by means of continuum- 
continuum transitions. Detailed studies of dynamic con- 
tributions to the breakup were explored in refs. [5U [53] 
and in several other publications which followed. The role 
of higher multipolarities (e.g., E2 contributions [51H5B] in 
the reaction ^Be(p, 7)^B) and the coupling to high-lying 
states has also to be investigated carefully. It has also 
been shown that the influence of giant resonance states 
is small [67]. 



D. Charge exchange reactions 

During core collapse, temperatures and densities are 
high enough to ensure that nuclear statistical equilib- 
rium is achieved. This means that for sufficiently low 
entropies, the matter composition is dominated by the 
nuclei with the highest binding energy for a given Y^.. 
Electron capture reduces Fg, driving the nuclear compo- 
sition to more neutron rich and heavier nuclei, including 
those with > 40, which dominate the matter compo- 
sition for densities larger than a few 10^" g cm~'^. As 
a consequence of the model applied in collapse simula- 
tions, electron capture on nuclei ceases at these densities 
and the capture is entirely due to free protons. To under- 
stand the whole process it is necessary to obtain Gamow- 
Teller matrix elements which are not accessible in beta- 
decay experiments. Many-body theoretical calculations 
are right now the only way to obtain the required matrix 
elements. This situation can be remedied experimentally 
by using charge-exchange reactions. Charge exchange re- 
actions induced in (p,n) reactions are often used to obtain 
values of Gamow- Teller matrix elements, B{GT), which 
cannot be extracted from beta-decay experiments. This 
approach relies on the similarity in spin-isospin space of 
charge-exchange reactions and /3-decay operators. As a 
result of this similarity, the cross section f7(p, n) at small 
momentum transfer q is closely proportional to B(GT) 
for strong transitions [55] , 



da 
dq 



{q = Q)=KND\Jar?B{a), 



(15) 



where if is a kinematical factor, is a distortion factor 
(accounting for initial and final state interactions). Jar 
is the Fourier transform of the effective nucleon-nucleon 
interaction, and B{a = F, GT) is the reduced transition 
probability for non-spin-flip, 

B{F) = {2J, + l)-\f\\Y,rlt\m\ 

k 
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and spin-flip, 

i3(GT) = (2J, + l)-i|(/||^a,r(±)||z)p, 

fc 

transitions. 

Eq. [15) valid for one-step processes, was proven to 
work rather well for (p,n) reactions (with a few excep- 
tions). For heavy ion reactions the formula might not 
work so well. This has been investigated in refs. [69H7T] . 
In Ref. [55] it was shown that multistep processes involv- 
ing the physical exchange of a proton and a neutron can 
still play an important role up to bombarding energies of 
100 MeV/nucleon. Refs. [70l [Tl] use the isospin terms 
of the effective interaction to show that deviations from 
the Taddeucci formula are common under many circum- 
stances. As shown in Ref. [72], for important GT tran- 
sitions whose strength are a small fraction of the sum 
rule the direct relationship between cr(p, n) and B{GT) 
values also fails to exist. Similar discrepancies have been 
observed |73j for reactions on some odd-A nuclei includ- 
ing ^'^C, "'^^N, '^^Cl, and '^^K and for charge-exchange 
induced by heavy ions [TT] [TJ. Undoubtedly, charge- 
exchange reactions such as (p,n), (^He,t) and heavy-ion 
reactions (A,A±1) can provide information on the B{F) 
and B[GT) values needed for astrophysical purposes [75] . 

E. Knock-out reactions 

Exotic nuclei are the raw materials for the synthesis 
of the heavier elements in the Universe, and are of con- 
siderable importance in nuclear astrophysics. Modern 
shell-model calculations are also now able to include the 
effects of residual interactions between pairs of nucleons, 
using forces that reproduce the measured masses, charge 
radii and low-lying excited states of a large number of 
nuclei. For very exotic nuclei the small additional stabil- 
ity that comes with the filling of a particular orbital can 
have profound effects upon their existence as bound sys- 
tems, their lifetimes and structures. Thus, verifications 
of the ordering, spacing and the occupancy of orbitals 
are essential in assessing how exotic nuclei evolve in the 
presence of large neutron or proton imbalance and our 
ability to predict these theoretically. Such spectroscopy 
of the states of individual nucleons in short-lived nuclei 
uses direct nuclear reactions. 

The early interest in knockout reactions came from 
studies of nuclear halo states, for which the narrow mo- 
mentum distributions of the core fragments in a qualita- 
tive way revealed the large spatial extension of the halo 
wave function. It was shown |76j that the longitudinal 
component of the momentum (taken along the beam or 
z direction) gave the most accurate information on the 
intrinsic properties of the halo and that it was insensitive 
to details of the collision and the size of the target. In 
contrast to this, the transverse distributions of the core 
are significantly broadened by diffractive effects and by 
Coulomb scattering. For experiments that observe the 



nucleon produced in elastic breakup, the transverse mo- 
mentum is entirely dominated by diffractive effects, as il- 
lustrated [77] by the angular distribution of the neutrons 
from the reaction ^Be(^^Be,^''Be-|-n)X. In this case, the 
width of the transverse momentum distribution reflects 
essentially the size of the target [78] . 

Most practical studies of medium corrections in 
nucleon-nucleon scattering are carried out by consider- 
ing the effective two-nucleon interaction in infinite nu- 
clear matter. This is known as the G-matrix method, 
an is obtained from a solution of the Bethe-Goldstone 
equation 

(k|G(P,pi,p2)|ko) = (k|v^Ar|ko) - 

f <fik' (k|vJVjv|kOQ(k^P,pl,p2)(k'|G(P,pl,p2)|ko) 
J (2^)3 ^;(P,k')-^o-^e 

(16) 

with ko, k, and k' the initial, final, and intermediate 
relative momenta of the NN pair, k = (ki — k2)/2 and 
P = (ki -|- k2)/2. If energy and momentum is conserved 
in the binary collision, P is conserved in magnitude and 
direction, and the magnitude of k is also conserved, ^mn 
is the nucleon-nucleon potential. E is the energy of the 
two-nuclcon system, and Eo is the same quantity on-shell. 
Thus £;(P, kj = e(P + k) + e(P - k), with e the single- 
particle energy in nuclear matter. It is also implicit in 
Eq. [16] that the final momenta k of the NN-pair also lie 
outside the range of occupied states. 



T 1 1 — I — I I I I I 1 1 — r 
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10 20 50 100 200 500 

E|ab (MeV/nucleon) 

FIG. 4: Total knockout cross sections for removing the / = 
halo neutron of ^^C, bound by 1.218 MeV, in the reaction 
^Be(i^C,"C9s). The solid curve is obtained with the use 
of free nucleon-nucleon cross sections. The dashed curve in- 
cludes the geometrical eflfects of Pauli blocking. The dashed- 
dotted curve is the result using the Brueckner theory, and the 
dotted curve is a phenomenological parametrization. 

In Ref. [89J the numerical calculations have been per- 
formed to account for the geometric effect of Pauli block- 
ing. A parametrization has been devised which fits the 
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numerical results. The parametrization reads 



\ Po J \ PPO J 



2.75 



1 - 



, if S > 46.27p2/3 



23tJ^' if £; < 46.27p2/3 

(17) 

where E is the laboratory energy in MeV, pi is the local 
density of nucleus i, /?< = niin(pi,/92) and p — {pi + 
P2)/po, with Po — 0.17 fm^^. 

The Brueckner method goes beyond a treatment of 
Pauli blocking, and has been presented in several works, 
e.g. in Ref. [90| I91j. where a simple parametrization was 
given, which we will from now on refer as Brueckner the- 
ory. It reads (the misprinted factor 0.0256 in Ref. [HI] 
has been corrected to 0.00256) 



250 MeV/nucleon 
5Be("Be,'"Be) 




-40 40 

p, (MeV/c) 



FIG. 5: Longitudinal momentum distribution for the residue 
in the ^Be("Be,^°Be), reaction at 250 MeV/nucleon. The 
dashed curve is the cross section calculated using the NN 
cross section from the Brueckner theory and the solid curve 
is obtained the free cross section. 



31.5 + 0.092 20.2- £;° 



1 + 0.0034£;i-5ip2 
1 + 21.55pi-34 



1.05 ^3 



= 23.5 + 0.00256 (18.2 -i?-)- ' + "^'^'fj^,/' 

(18) 



A modification of the above parametrization was done 
in Ref. |92j . which consisted in combining the free nu- 
cleon nucleon cross sections parametrized in Ref. |93j 
with the Brueckner theory results of Ref. [501 IHI] • 

To test the influence of the medium effects in nucleon 
knockout reactions, we consider the removal of the I = 
halo neutron of ^^C, bound by 1.218 MeV, and the Z = 
neutron knockout from "^^Ar, bound by 17.06 MeV. The 
reaction studied is ^Be(^^C,^'*Cgs). The total cross sec- 
tions as a function of the bombarding energy are shown 



in figures IV E The solid curve is obtained with the use 
of free nucleon- nucleon cross sections. The dashed curve 
includes the geometrical effects of Pauli blocking. The 
dashed-dotted curve is the result using the Brueckner 
theory, and the dotted curve is the phenomenological 
parametrization of the free cross section. 

In figure [5] we plot the longitudinal momentum 
distributions for the reaction ^Be(^^Be,^'^Be), at 250 
MeV/nucleon [89 . The dashed curve is the cross section 
calculated using the NN cross section from the Brueck- 
ner theory and the solid curve is obtained the free cross 
section. One sees that the momentum distributions are 
reduced by 10%, about the same as the total cross sec- 
tions, but the shape remains basically unaltered. If one 
rescales the dashed curve to match the solid one, the dif- 
ferences in the width are not visible [Ml. 



V. CONCLUSIONS 

There were many questions not addressed in this re- 
view, such as the role of central nucleus-nucleus colli- 
sions in determining phase transition, equation of state, 
and a quark-gluon plasma, all topics or relevance in as- 
trophysics. The review was more focused on the role of 
short-lived, exotic nuclei. The important scientific ques- 
tions to be addressed both experimentally and theoret- 
ically in nuclear physics of exotic nuclei with relevance 
for astrophysics comprise: (a) How do loosely-bound sys- 
tems survive and what are the general laws of their for- 
mation and destruction? (b) Are new types of radioac- 
tivity possible? (c) Are new types of nuclear symmetry 
and spatial arrangement possible? (d) What are the lim- 
its of nuclear existence? (e) How do the properties of 
nuclear matter change as a function of density, temper- 
ature and proton-to-neutron ratio? (f) How do thermal 
and quantum phase transitions occur in small systems? 
(g) What determines the shape and symmetry properties 
of an exotic nucleus? (h) How does quantum tunneling 
of composite particles occur in the process of reactions 
and decay? (i) What are the manifestations of funda- 
mental forces and symmetries in unusual conditions? (j) 
How were the elements heavier than iron formed in stellar 
explosions? (k) How do rare isotopes shape stellar explo- 
sions? (1) What is the role of rare isotopes in neutron 
stars? These questions provide extreme challenges for 
experiments and theory. On the experimental side, pro- 
ducing the beams of radioactive nuclei needed to address 
the scientific questions has been an enormous challenge. 
Pioneering experiments have established the techniques 
and present-generation facilities have produced first ex- 
citing science results, but the field is still at the beginning 
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of an era of discovery and exploration that will be fully 
underway once the range of next- generation facilities 
becomes operational. The theoretical challenges relate 
to wide variations in nuclear composition and rearrange- 
ments of the bound and continuum structure, sometimes 
involving near-degeneracy of the bound and continuum 
states. The extraction of reliable information from ex- 
periments requires a solid understanding of the reaction 
process, in addition to the structure of the nucleus. In as- 



trophysics, new observations, for example the expected 
onset of data on stellar abundances, will require rare- 
isotope science for their interpretation. 
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